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Abstract 

We use an effective field theory to compute low-energy nucleon- 
deuteron scattering. We obtain the quartet scattering length using low 
energy constants entirely determined from low-energy nucleon-nucleon 
scattering. We find a t h = 6.33 fm, to be compared to a exp = 6.35±0.02 
fm. 
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There has been considerable interest lately in a description of nuclear 
forces from the low-energy effective field theory (EFT) of QCD. (For a review, 
see Ref. [[]]].) Following a program suggested by Weinberg 0, the leading 
components of the nuclear potential have been derived || and a reasonable 
fit to two-nucleon properties has been achieved [|J. The correct formulation 
of the nuclear force problem within the EFT method is important because 
it will allow a systematic calculation of nuclear properties consistently with 
QCD. One would like, for example, to devise a theory of nuclear matter 
rooted in a hadronic theory that treats chiral symmetry correctly and yields 
the well-known few-nucleon phenomenology. One hopes that after a number 
of parameters of the EFT are either calculated from first principles or fitted 
to a set of few-nucleon data, the theory can be used to predict other reactions 
involving light nuclei and features of heavier nuclei. 

However, some issues concerning renormalization in this non-perturbative 
context and fine-tuning in the two-nucleon S"-waves have been raised in Refs. 
[|, U| and are still not fully understood [[/J. The fine-tuning necessary to 
bring a (real or virtual) bound state very close to threshold generates a scat- 
tering length a much larger than other scales in the problem. At momenta of 
0(1 /a), mesons can be integrated out and the characteristic mass scale [i of 
the underlying theory controls the size of the other effective range parameters; 
for example, the effective range ro ~ 2//x. Once the leading order contribu- 
tions, which give rise to a, are included to all orders, the EFT at momenta 
0(1/ a) becomes an expansion in powers of l/(ap). Kaplan || has noticed 
that the interactions that generate a non-zero r can also be resummed by 
the introduction of a baryon number two state of mass A = 2/Ma\r \, which 
in lowest order in a derivative expansion couples to two nucleons with a 
strength g 2 /^ = l/M 2 |ro|. In Ref. || it was shown how this works in the 
two-nucleon 1 S'o channel. Analogous considerations hold for the 3 Si channel, 
where they are similar to the old quasi-particle approach of Weinberg || . 

In this paper we consider the application of these ideas to the three- 
nucleon system. Our goal here is to calculate some of the three-nucleon pa- 
rameters that are dominated by the leading interactions in the EFT without 
pions. We show in particular that the quartet scattering length in neutron- 
deuteron scattering can be predicted once the EFT is constrained by low- 
energy two-nucleon data. Such an attempt to a model-independent or "uni- 
versal" approach is not a new idea; it permeates for example the work of 



Efimov (see, e.g., Hj) and Amado (see, e.g., pi). However, as we will show, 
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the EFT formulation is much easier to implement, from both conceptual and 
practical standpoints. 

For momenta of order 1/a (the momentum scale relevant for zero-energy 
Nd scattering), we can integrate out mesons and consider an EFT with 
only nucleons N. Interactions are then described by a tower of nucleon 
contact operators with an increasing number of derivatives. Amplitudes in 
leading order are given by a zero-range four-nucleon interaction iterated to 
all orders. Corrections come in powers of l/(am n ). The next two orders in 
this expansion, l/(am n ) ~ r^jla and l/(am n ) 2 ~ (ro/2a) 2 , stem from one 
and two insertions of a two-derivative four-nucleon operator giving rise to 
a non-zero r®. It is advantageous to sum all the contributions coming from 
this operator, which can be easily done because they appear in a geometric 
series. The resulting interaction is equivalent to the s-channel propagation of 
a dibaryon, and therefore can be obtained more directly by the introduction of 
a dibaryon field. Assuming naturalness, only higher orders depend on further 
NN scattering parameters — such as the shape parameter, which contributes 
at 0(l/(am 7r ) 3 ) — and three-nucleon forces, which start at 0(l/(am n ) 4 ). 

Since in both 1 = and 1=1 S-wave two-nucleon channels we observe 
(one real, one virtual) bound states near threshold, we consider two dibaryon 
fields, T (D) of spin zero (one) and isospin one (zero). The most general 
Lagrangian invariant under parity, time-reversal, and small Lorentz boosts 
is 

C = NHid + ^ + ...)N 

+Tt.(-id -^ + A T + ...)T + &.(-id -^ + A D + ...)D 
_^I( T t . Na 2 rr 2 N + h.c.) - ^(Et ■ Nr 2 aa 2 N + h.c.) + . . . (1) 

Here the A T>D and g T D are undetermined parameters and ". . ." stands 
for higher order terms. (Note that the effects of non-derivative and two- 
derivative four-nucleon terms can be absorbed into a redefinition of At,d 
and gT,D and higher order four-nucleon terms.) 

In this non-relativistic theory all particles propagate forward in time, 
nucleon tadpoles vanish and, as a consequence, there is no dressing of the 
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Figure 1: Dressed dibaryon propagator. 



nucleoli propagator, which is simply 



Sn(p) 



i 



(2) 



The propagators for dibaryons are more complicated, because of the coupling 
to two-nucleon states. The dressed propagators consist of the bubble sum 
in Fig. [I], which amounts to a self-energy contribution proportional to the 
bubble integral. This integral is proportional to the (large) mass M, and it 
is this enhancement that gives rise to non-perturbative phenomena and leads 
eventually to the existence of bound states. The integral is also ultraviolet 
divergent and requires regularization. Introducing a cut-off A we find a lin- 
ear divergence oc A, a cut-off independent piece which is non-analytic in the 
energy, a term that goes as A -1 and terms that are higher order in A -1 . The 
first and third terms can be absorbed in renormalization of the parameters 
of the Lagrangian ([]]); in what follows we omit a label R that should be 
attached to these parameters, i.e., At,d and gr,D stand for the renormal- 
ized parameters. Higher order terms are neglected because they are of the 
same order as interactions in the ". . ." of the Lagrangian ([!]). A dibaryon 
propagator has therefore the form 



Note that such a dressed propagator has two poles at p° = p 2 /AM — B, p° = 
p 2 /AM — B deep and a cut along the positive real axis starting at p° = p 2 /AM. 

The AW amplitude can now be obtained directly from Sn(p) as in Fig. 
0. In the center-of-mass, the on-shell / = 0, J = 1 S-wave amplitude at an 
energy E = k 2 /M is 
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Figure 2: NN amplitude. 



which is exactly equivalent to the effective range expansion. An analogous 
result holds for the 1=1 S-wave. The four parameters At,d and gr,D can 
then be fixed from the experimentally known scattering lengths and effective 
ranges. The NN amplitude has shallow poles at B ~ 1/Ma 2 which are 
associated with the deuteron in the 3 Si channel and with the virtual bound 
state in the 1 Sq channel. The effective theory has also an additional deep 
bound state in each channel at B deep ~ 4/ Mr 2 ., which is outside the range of 
validity of the EFT. 

From the triplet parameters 3 a = 5.42 fm and 3 r = 1.75 fm ID] we find 



A/) = 8.7 MeV and g 2 D = 1.6 ■ 10 -3 MeV. The resulting deuteron binding 
energy is B = 2.28 MeV. From the singlet parameters l a pp = —17.3 fm, 
1 anp = —23.75 fm, 1 a nn = —18.8 fm, 1 r 0pp = 2.85 fm, 1 r Qnp = 2.75 fm, and 
Von,, = 2.75 fm fljl we find the averages A T = —1.5 MeV and g\ = 1.0- 1CT 3 
MeV. 

With the parameters so determined, we turn now to possible predictions 
in low-energy nucleon-deuteron scattering. For simplicity we restrict our- 
selves to scattering below the deuteron break-up threshold, where the S- 
wave is dominant. There are two S-wave channels, corresponding to total 
spin J = 3/2 and J = 1/2. In the quartet only D contributes while in the 
doublet T also appears. The Nd scattering amplitude T/vd from the same 
interactions is given by the diagrams in Fig. ^|, which can be summed up 
by solving an integral equation in the quartet and a pair of coupled inte- 
gral equations in the doublet channel. The diagrams under consideration are 
power-counting finite, but this does not preclude the existence of relevant 
contact interactions between the nucleon and the dibaryons. Pion exchange 
that would generate such interactions can be expected to be larger for the 
1 = 1 dibaryon T, and therefore predominantly affect the J = 1/2 chan- 
nel. We will return to the doublet case in a future publication. Here we 
study the quartet channel expected to be much less sensitive to the details 
of the physics of momenta of 0(m n ), since the wave function on this channel 
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Figure 3: Nd amplitude. 

vanishes by symmetry when the three particles are at the same point. 

An enormous simplification comes about because the s-channel interac- 
tion due to the dibaryon is both local and separable. This allows us to write 
a simple integral equation that sums all the graphs in Fig. [| Performing the 
integration over the time-component of the loop 4-momentum, we find that 
the conveniently normalized on-shell amplitude as a function of the initial 
(final) center-of-mass 3-momentum k (p) satisfies 



3(P " 2 k2 K^J\{p-P)+MB-^MB) 



8M 2 g 2 



D 



t(p,k) 



p 2 — k 2 — ie 



(5) 



dH 



1 



t(l,k) 



(p + k/2) 2 + MB J (2n) 3 [2 + i . p + p - ^k 2 + MB I 2 - k 2 



Note that all terms in a perturbative expansion of t in (|5|) are of the same 
order (~ l/yMB). It is straightforward but tedious to show that the wave 



function (27r) 3 5(p — k)+t(p, k)/(p 2 — k 2 



■ie 



corresponding to a scattering so- 



lution indeed satisfies the Schrodinger equation derived from the Lagrangian 
©■ 

At zero energy {k — > 0) only the S-wave, depending on the magnitudes 
of momenta, contributes to the scattering, and we can perform the angular 
integration directly. It is also convenient to normalize all quantities to V MB. 
Defining x = p/V MB , 



a(x) 



4tt ^y/MB 



,0), 
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Figure 4: Function a(x) for 77 = 0.40 without (solid line) and with (dashed 
line) cut-off. 



and introducing 



Eq. 



becomes 



< f , ,x 2 + z 2 + 1 + xz. 

F{x,z) = —In , 

xz x z + z^ + 1 — xz 



(7) 



—r] + 



1 + ^/1 + fx 2 
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dzF(x, z)a(z) 



Note that there is only one parameter 77 = 2tt a/ M73 / M 2 g 2 D = 3 r \/ MB/2 = 
0.40 in this equation. The value of the function a(x) at a; = gives the Nd 
scattering length in units of 1/ V MB . The same equation was previously 
obtained and solved in the zero-range limit (rj — > 0) [|i~3|j . 

We have solved Eq. (|8|) numerically for rj = 0.40 by the Nystrom method 
[|i~4|1 . The solution a(x) is plotted as the solid line in Fig. ^. The pole in 
a(x) around x ~ 4.4 is associated with the spurious deep two-body pole. Its 
presence allows intermediate states where two nucleons fall into this deep 
state while the other has extra energy. This means that the outgoing wave 
has an additional component, a pole at the momentum corresponding to this 
additional process. The interesting point is that even though the effective 
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theory makes nonsensical predictions outside its domain of validity, like the 
existence of this new state in Nd scattering, the low-x part of the curve 
is insensitive to the large-x behavior, and the prediction for the scattering 
length is sensible. In order to demonstrate this more explicitly we have also 
solved Eq. (H) with a cut-off two-nucleon amplitude without the deep pole. 
For a cut-off of 150 MeV we obtain the broken line in Fig. |j. 

The quartet scattering length is 4 a = -a(0)/vMB. For rj = (and 
B fixed), we reproduce the result 4 a = 5.09 fm of Ref. |I5fl. Taking into 



account the finite range (77 = 0.40) we obtain (Fig. |j) 4 a = 6.33 fm with 
an uncertainty from higher orders of ~ ±0.10 fm. This result obtained with 
no free parameters is in very good agreement with the experimental value of 
4 a = 6.35 ±0.02 fm WB. 
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